Pion stability in a hot dense media 
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Pions may remain stable under certain conditions in a dense media at zero temperature in the 
normal phase (non pion superfluid state). The stability condition is achieved when the in- media pion 
width vanishes. However, thermal fluctuations will change this stable regime. For low temperature 
pions will remain in a metastable state. Here we discuss the different possible scenarios for leptonic 
pion decays at finite temperature, taking into account all the different chemical potentials involved. 
The neutrino emission due to pions in a hot-dense media is calculated, as well as the coolig rate of 
a pion-lepton gas. 
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The study of pions in a superfluid state has called 
the attention of physicist during many years in different 
frameworks |ll4lOl|. In particular, compact stars may pro- 
vide a natural scenario for such a state of matter. How- 
ever, in the absence of a stable superfluid state, pions in 
the normal phase could behave as metastable particles. 
In fact, for high densities, an appropriate combination 
of chemical potentials can avoid the normal pion decay 
process. Indeed, if in cold matter the leptons states are 
completely filled up to the Fermi level, there will be no 
possible allowed final states for the decay of charged pi- 
ons, unless thermal fluctuations are present. Pion decay 
properties in a hot and/or dense media have been vastly 
studied pH - flq . The present work is focused on the ex- 
istence of possible pion metastable states in the normal 
phase at high densities. As a consequence, we will study 
how these metastable states affect the neutrino emission 
rate. The cooling rate for a pion-lepton gas is also deter- 
mined. The temperature considered is much less than the 
pion mass. We will consider baryon chemical potential 
values lower than ^ 1 GcV, in order to avoid other pos- 
sible media effects like color superconductivity. The pos- 
sible scenario for such a process are protoneutron stars 

El, El. 



I. LOW ENERGY QCD AT FINITE DENSITY 

As the 99.9877% of the charged pions decay into muons 
and muonic neutrinos, we will refer speciflcally to this 
process. Our discussion will be based on the following 
low-energy Lagrangian: C = Cf^ 



£p^, where 



E ^/ 



- M/70 



V, 



t/-/ 



(1) 



/=M,''m 



corresponds to the lepton free Lagrangian, fif being the 
associated lepton chemical potential. The second term 
in the low-energy Lag rangian corresponds to the 0(p)^ 
chiral Lagrangian [20|, [21 1 



C. 



G 



The F and G terms are the tree level pion decay con- 
stant and the tree level chiral condensate, respectively. 
The U fields contain the pion fields as C/ = exp(i7r"r°//) 
and the covariant derivative D^JJ — d^U — ir^U + iUl^ 
includes external right and left currents. When M. = 
diag(?7i„, TOd) we break explicitly the chiral symmetry. 
In this article we are not interested in computing pion 
radiative corrections. Therefore, higher order chiral La- 
grangian terms will not be considered. However, those 
corrections have been already calculated at finite tem- 
perature and isospin chemical potential [16i|, and there- 
fore, if we want to incorporate these contributions, it is 
enough to replace the masses, decay constant and chem- 
ical potential terms by temperature and isospin chemical 
potential dependent dressed terms 22 1. 

We will use the effective Fermi model for the lep- 
tonic weak coupling. The leptonic weak currents and 
the isospin chemical potential are introduced by setting 
the external currents in the chiral Lagrangian as 



la = GfH'u^Ic 

1 



(l-75)-0^r +^^7a(l -75)V'i/„ 



+ 7:A*/T 5aG 



(3) 



(4) 



with Hi = fJ.u ^ fJ-d is the isospin chemical potential and 
where we use the following combination of Pauli matri- 



ces: r 



-7^(r^ ± ir^). We will concentrate only on 
I A*/ 1 < "^TT- In the superfluid 



71^ 



tr{D^UyDU'' + -tr{U''M + M^U). (2) 



the normal phase where |/z/| < m 

phase, \fii\ > mTT, one of the charged pions condenses, 

and therefore, another treatment is needed @, Q- 



Baryon chemical potential 

Two different regions in the literature have been con- 
sidered when the Baryon chemical potential /is = | (/i^ -I- 
fid) is introduced in the frame of chiral perturbation the- 
ory: 

Small fiB ■ The baryon chemical potential can be taken 
as 0{p) in the power counting in chiral perturbation the- 
ory. In this case /is appears in the Wess-Zumino-Witten 



anomalous term, which turns out to be relevant only in 
higher order radiative corrections |23| . 

Very high fis- In this case an expansion in powers 
of fiB~^ is performed (asymptotically infinite chemical 
potential). Here, effects of color superconductivity and 
color flavor locked phase are present, due to the ap- 
pearance of diquark pairing [1^, [24| - |26| . We can con- 
struct an effective Lagrangian including hb in the ab- 
sence of diquark effects, by considering effective F{ij,b) 
and G{^b) constants in the chiral Lagrangian in Eq. ^ 
and introducing a Lorentz symmetry breaking term: 
{Df.UYDUt' -j> {DqUYDoU - v^{DUY ■ DU. Results 
obtained in the frame of the Nambu-Jona-Lasinio model 
show that for ^j-b ^^ GeV the -FjG and v parameters do 
not suffer significant changes [3, [ij, IT3] ■ We can neglect 
then baryon chemical potential effects in this work. 

Once the chiral Lagrangian has been expanded in 
terms of the pion fields, the canonical quantization proce- 
dure is the standard one, keeping in mind that the energy 
of the charged pions and lepton fields are shifted due to 
the chemical potentials. Then, we proceed to calculate 
the decay width of the charged pions and the correspond- 
ing neutrino emissivity. 



II. CHARGED PIONS DECAY WIDTHS 

The decay width for charged pions, including finite 
temperature and density effects, is given by 

r.± = ^ y dq^± dk^^ (27r)4^(4) (p _ ^ _ fc) 

x\M±\''[l-nFiqo)-nFiko)], (5) 

where the on-shcU pions are in the rest frame, p ~ (m^ =p 
/i/,0), and where np{z) — [e^l^ -I- f)~^ is Fermi-Dirac 
distribution. The phase space measure is defined as 



dkp = 



_(£k_ 



e{k„+^lp)6{{ko+^ipf-k^-ml), (6) 



where P stands for the different particles involved: pions, 
fermions and antifermions. The transition probability 
matrix was abbreviated as 



lA^iP^EK'^^'^;^!^'"'!^^)! 



spin 



The corresponding chemical potentials are /i^± = ifJ-i, 
and iifT — ifJ-f where /^ (/^) denotes a fermion (an- 
tifermion) . 

This definition of the decay width corresponds to the 
imaginary part of the thermal one-loop weak interaction 
corrections to the pion propagator. The decay width in- 
cludes the decay of pions into leptons as well as recombi- 
nation of leptons and neutrinos. By considering massless 
neutrinos, the decay width for charged pions is then 



r^± = Ftt 



f -n_F(e^±) -n_F(e^q:) 0{m^ ^ S^i ~ m^_,) 



with 



1 — rai/rn^ 



-n Jtt'-'F 2 fi 2/21 
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2(m„ T Sji) 
2(m^ =F Sfi) 
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(8) 

(9) 
(10) 

(11) 
(12) 



Fjr being the vacuum pion decay width. 

For the non chemical equilibrium case, when \Sfi\ > 
mTr—Tn^, one of the charged pions remains stable since its 
decay width vanishes, as can be seen from the Heaviside 
function in Eq. (|S]). We will consider then, in the non- 
equilibrium case, that \6^\ < to^ — m^ « 34 MeV. Here 
we used m^ = f 39.6 MeV and m^ = f05.6 MeV. 

Stable states at zero temperature, i.e. those where 
their decay width vanish, will now develop a small decay 
width due to thermal fluctuations. From Eq. ([5]), con- 
sidering that nF{x) — )■ 9{—x) when T -> 0, we find two 
metastable cases: 



met ast able tt 



metastable 



Mm > TTTZ — ):.A : (13) 



Mi^^, > 



2(?7l7r + SfJ.) 

(m^ - (5^)2 - ml 
2(m7r — Sfi) 



(14) 



Leptonic and beta equilibrium 

In terms of quark degrees of freedom, the beta equi- 
librium condition is Hd = IJ'u + IJ-e — y^v^ ■, or, in terms 
of the isospin chemical potential, ^i + fj,e ^ /^jy^ =0. If 
we consider a degenerate gas, which is the case for lep- 
tons in compact stars, the condition /i^ — /i;^^ = Me — Mi^e 
arises in order to equilibrate the Fermi levels [1^. As 
a consequence of this leptonic chemical equilibrium, the 
beta-equilibrium condition will produce S^ = 0. 

From Eqs. (|13p and (fH|) , we can see that metastable 
states in beta equilibrium will occur for jjl^ > /i* for 7r~ 



(7) mesons and //^^ > fj,* for 7r+ mesons, where 



m;: 



2TO7r 



2TO7r 



109.74 MeV, 



29.9 MeV. 



(15) 
(16) 



Note that in beta equilibrium, the tt^ meson condenses 
if the lepton chemical potential is high enough such that 
fJ-fj. ~ fJ-iy > 'm-TT- On the other side, if the neutrino chem- 
ical potential is high enough, such that /i;y, — /^^ > m-^, 
then a 7r+ meson condenses [42]. 

In order to explore the phenomenological consequences 
of these pion metastable states, we will discuss next the 



neutrino emission and cooling rate of a pion-lepton gas 
in leptonic- and beta- equilibrium. Our discussion will 
show in a clear way the influence of metastable states on 
the cooling rate. 



III. NEUTRINO EMISSIVITY 

Neutrino emission is perhaps the most relevant phe- 
nomena associated to the temperature evolution of com- 
pact stars. The neutrino emissivity e is defined as the 
energy loss through neutrino emission per unit time and 
unit volume. For the decaying pions, the neutrino emis- 
sivity must include the probability of finding a pion in 
the media as well as the Pauli blocking for the emerging 
muons 

e,± = J dp^^dq^^dk^i {2n)^S^^\p -q-k) 

x\M±\^ ko nB{po)[l - npiqa)], (17) 

where nB{z) = [e^l^ — 1)^^ is the Bosc-Einstein distri- 
bution and the other terms were defined in the previous 
section. 

Hereafter, we will consider the degenerate case where 
M/i ^ '^/i s-iid ^Iso beta-equilibrium. As we discussed in 
the last section, these two assumptions imply J/i = 0. 

The neutrino emissivity, then, is given by 



V^ml 



^<^ Tr'^{m'^ — m?) 



rK 
dE^^ / dEfj_ (E'tt — -E^, 



xnB(^^T(AV -Mi^J)[l -'^f(£'mTMm)]' (18) 
where the limits for the muon energy integral are 



E^. 



1 



2111,1 



{ml 



)E. ± {ml 



\)VeT 



mt, 



(19) 
The temperature and the muon chemical potential 
tend to favor the anti-neutrino emission due to the 7r~ 
decay. The ns factor in Eq. (|18p , which gives the prob- 
ability of finding a n^ meson, grows as function of tem- 
perature and lepton chemical potential. On the other 
hand, the (1 — np) factor gives the probability of finding 
an accessible state for the emerging muon. Only thermal 
fluctuations will conspire against the Pauli blocking, al- 
lowing then the decay of 7r~ . The probability of finding a 
7r+ meson becomes smaller for higher chemical potential, 
being suppressed by the Bose factor. 

We are interested to extract the leading terms for the 
low temperature behavior of the emissivity. Since we 
arc considering a degenerate gas (/i^ > rn.^), the main 
contribution is given by the emission of antineutrinos, 
being the neutrino emission highly suppressed by an ex- 
ponential factor as we mentioned previously. In order 
to extract the main contribution to the total emissivity, 
we expand Eq. (J18p in the low temperature region at 
the leading order. The Fermi-Dirac distribution becomes 




FIG. 1. Efj^ as a function of Ett, indicating three different 
regions for /i^ (horizontal lines) which determine the inte- 
gration regions of the Et^ variables for the emissivity at low 



temperature, ruled by the condition E^^ > fi^ in Eq. pCT 



then nF{Ef, - ^^) « ^(/i^ - Efj_), obtaining 
e,^ « / dE^[g+ e{E+ - fi^) + g^ 0{E- - /.^)], (20) 



with 



9± 



± 



T^mlnB{E^ - Hf, + H^J 
7r2(m2 -m2) 



EttE,, 



2^^ 



(21) 

We can separate the antineutrino emissivity in three dif- 
ferent regions, depending on the value of the lepton chem- 
ical potential. Fig. [1] shows the functions E^ plotted as 
a function of Et^. On the vertical axes, the three regions 
are indicated for specific values of the lepton chemical 
potential. The integration limits in equation (|20p are 
determined by the condition E^ > /i^. 

In region I, where /i^ < m^, the argument of the theta 
functions is always positive, then 



/•OO /"OO 

ep^ (I) = / dEyrg+ + / dE^g^ 



(22) 



In region II, where to^ < A'p < t^-ln the function E'^ is 
positive for all values of E.^. However, from Eq. ((20|) . the 
condition E~ > fi^ will exclude a region in the integral: 

/>00 />771_ />00 

ep^(II)=/ dE^g++ dE^g^+ dE^g^,{23) 



where 



m± = 



2ml 



{ml + m^)/!^ ± {ml - to^)^, 



^^l - ml 



are the solutions of the equation E^{Et^) 



(24) 
fJ-ii, giving 



as a result _E^ = m±. This is an intermediate region 
between low and high lepton chemical potential. 

In region III, where /i^ > f^*ai the condition Ef^ > /i^ 
will exclude some values in the integration limits; 



/•CXD /"OO 

ep^(III)= / dE^g++ / dE^g^ 

Jm- Jm^ 



(25) 



where m± was previously defined in Eq. ([24)) and cor- 
responds to the solutions of the equation E^{Et^) — /i^, 
giving E+{m^) = E-{m+) = /i^. 

As we can see, the integrals above can be written in 
the form 



dE^f{Ej,)nBiEj, - fi^ + fi^) 



(26) 



where the integrand / ns — g± in Eq. (|2ip. and where 
m stands for 771^1, m±. In order to extract the leading 
terms in the low temperature region, ii m > fi^ — fj,,y , we 
can expand the Bose-Einstein distribution, and through 
an appropriate change of variables, we find 

00 
I=Y^ e-'^"(™-'^^+'''')T"F„,„(T, m) (27) 



n=l 



e-'^('"-'^"+'''')r"F„,i(0,m) 



with 



Fa,niT,m) 






(28) 



(29) 



and with a such that the last integral remains finite in the 
limit T ^f Q. Due to the exponential factor, the integrand 
in the above equation will be dominated by low a;-values. 
If the condition to — /i^, — /Xj^ > T is satisfied, we can keep 
only the first term n = 1 in the series. If ttt, « /i^ — /i^ ^ , 
we need to sum the whole series. 

As a result, the low temperature behavior of the emis- 
sivity becomes 



A T'^I'^ Q-{rn,-fj.^ 



BT e 



-(m_-/j^+p, 






for 



Mm 



for /x^ > ^* 



where 



A = r^TO^(l - ml/ml){2T:m^) 



B = r^TO^ ^ 

7r(TO^ 



-(27r?7l7r 



-3/2 
-1 



(30) 

(31) 

(32) 



The intermediate region m^ < M/i < /^m will be a com- 
bination of terms '^ T^/"^ and ^ T. When the chemical 
potential grows, the linear term in T starts to dominate. 
The neutrino emissivity, at low temperature, is 
strongly suppressed by a term exp[— (TO,^-|-/i^ — ^^^)/T]. 
We are interested in the high chemical potential region 
since it increases the emissivity in the low temperature re- 
gion. Our numerical analysis suggests that these approx- 
imations will be valid for temperatures less than 50 MeV. 



IV. COOLING RATE 

In order to explore some effects in the metastable re- 
gion of the 7r~, Eq. (IT^ . we will now calculate the 
cooling rate due to muonic neutrino emission for a pion- 
lepton gas, at constant volume and charge. We will 



consider that /i^^ =0, which means that all the neu- 
trinos will escape from the gas. We also consider /3- 
equilibrium, lepton-equilibrium and neglect the process 
of neutrino emission through muon decay. In other words 
—fij — fifi — iJ-e- This model is a simplification, eventu- 
ally valid as isolated bubbles inside the nuclear media of 
compact stars, although finite volume effects should be 
taken into account. 

The cooling time t is defined as 



t^tn 



'-^dT, 



(33) 



To 



where to and Tq are the initial time and temperature, 
respectively, cy and e correspond to the specific heat 
and the emissivity, respectively. 

The specific heat per unit of volume is given by 



cv 



T fdS 



V \dT 



(34) 



V,N 



where S is the entropy, and N, in our case, is the 
charge number. In a non-interacting degenerated gas, at 
low temperature and high chemical potential, the charge 
number is dominated by fermions, and the leading term 
in the temperature expansion is constant: 



N : 



V 
3^ 



{^,l-mlf/^ + {^,l-mlf/^ 



(35) 



So, for such low temperature approximation, the constant 
N condition is equivalent to consider a constant /i. The 
specific heat per unit of volume can be written as 



Cv 






(36) 



where p is the pressure. 

The contribution to the specific heat per unit of volume 
for noninteracting fermions is 



CVf 



_9_ 

rp2 



(2^ 



(E - p.fnF{E ~ fi)nFin - E), (37) 



with g = 2. The formula for the specific heat per unit 
of volume, for bosons, is the same as for fermions but 
with g = —1 and changing np by ub- The relevant 
contribution comes essentially from electrons and muons, 
where the specific heat per unit of volume for a degen- 
erated fcrmion gas of mass m and chemical potential /i 
is cvj = ^ny' jj? — rn^T ii ^ > m. In the case of muons, 
when fie ~ mfj_, their contribution will be « 0.4(TO^r)^/^. 
Fig. [2] shows the temperature dependence, as func- 
tion of the logarithm of time in seconds, for three differ- 
ent values of the electron chemical potential: fie — tt^ii 
(muon degeneracy), fie — ^Ti-n (pion condensation) and 
another value in between: fi = 122.6 MeV. It can be 
seen from Fig. [5] that starting from an initial temper- 
ature, T = 5 MeV, the needed time to reach 1 MeV is 
extremely short, fraction of seconds, in the metastable 



> 



Me = 105.6 McV 
/Je = 122.6 McV 
;ie = 139.6 McV 





Log(t/ 



FIG. 2. Temperature as function of the logarithm of time in 
seconds, of a pion-lepton gas in /3 -equilibrium with zero neu- 
trino chemical potential. The initial temperature is 5 MeV. 



as well in the condensed case. On the other hand, it 
will take thousands of years to diminish the temperature 
from 1 MeV to fractions of MeV. The cooling process is 
not so fast if we consider values of the electron chemical 
potential lower than 109MeV. Notice that, in spite of the 
fact that /ig = 122.6 MeV is an average tween /ie = "t-^ 
and /ig ~ 'ni-n, the corresponding cooling time curve is 
notoriously closer to the beginning of the pion superfluid 
phase. 



V. CONCLUSIONS 

In this paper we have calculated the charged pions de- 
cay widths in dense matter at finite temperature, analyz- 
ing the pion metastable condition. We have calculated 
also the neutrino emissivity through the leptonic pion 
decay in /3-equilibrium for a degenerated system. We 
obtained the main contributions in the low temperature 
region for pions decaying into muons and niuonic neu- 
trinos. Finally, we estimated the cooling rate of a pion 
lepton gas in /3-equilibrium for three different values of 
the electric chemical potential. 

From our results, we argue that it is possible to find 
stable pions, even if they are not condensed, if the lepton 
chemical potential reaches a value higher than the muon 
mass, which is the case in compact stars. Under such 
conditions, the pions might only decay through thermal 
fiuctuations. The muonic neutrino emissivity will grow 
with the lepton chemical potential, varying from ~ T^' ^ 
to ^ T, both with an exponential suppressing term. The 
contribution to the cooling process in a neutron star, due 
to the emission of muonic neutrinos from pions, has not 
been much considered yet in the normal phase. In fact, 
from our estimation of the cooling time of the lepton 
pion gas, it can be seen that this process is relevant for 
temperatures higher than 1 MeV which corresponds to 
the cooling process of a protoneutron star. The pion 
superfluid case will be explored elsewhere. 
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